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ABSTRACT 



Aims. Guided by the duality of turbulence (random versus coherent motions), we seek coherent structures in the turbulent velocity 
field of molecular clouds, anticipating their importance in cloud evolution. 

Methods. We analyse a large map (40' by 20') obtained with the HERA multibeam receiver (IRAM-30m telescope) in a high latitude 
cloud of the Polaris Flare at an unprecedented spatial (11") and spectral (0.05km s~^) resolutions in the ^^C0(2 - 1) line. 
Results. We find that two parsec-scale components of velocities diff'ering by ~ 2 kms~\ share a narrow interface (< 0.15 pc) that 
appears as an elongated structure of intense velocity- shear, ~ 15 to 30 kms~^ pc~^ The locus of the extrema of line-centroid- velocity 
increments (E-CVI) in that field follows this intense-shear structure as well as that of the ^^C0(2- 1) high- velocity line wings. The tiny 
spatial overlap in projection of the two parsec-scale components implies that they are sheets of CO emission and that discontinuities 
in the gas properties (CO enrichment and/or increase of gas density) occur at the position of the intense velocity shear. 
Conclusions. These results disclose spatial and kinematic coherence between scales as small as 0.03 pc and parsec scales. They 
confirm that the departure from Gaussianity of the probability density functions of E-CVIs is a powerful statistical tracer of the 
intermittency of turbulence. They disclose a link between large scale turbulence, its intermittent dissipation rate and low-mass dense 
core formation. 

Key words. ISM: clouds, ISM: magnetic fields, ISM: kinematics and dynamics, turbulence 



1. Introduction 

Because it is supersonic, magnetized and develops in a multi- 
phase medium, interstellar turbulence is expected to differ from 
turbulence in laboratory flow experiments or in state-of-the- 
art numerical simulations, e.g. Chan al et a l. (2000) for rec ent 
experiments in gase ous helium and iMininni et alJ (l2006ah or 
lAlexakis et alJ (l2007l) for MHD simulations. Nonetheless, it may 
carry some universal properties of t urbulence, such as sp ace- 
time intermittency (for a review see lAnselmet et all 1200 ll) . Of 
particular interest to star formation, is the behavior of turbu- 
lence dissipation. In a series of papers, we have shown that 
the ^^C0(1 - 0) line-centroid- velocity increments (CVI) in 
translucent molecular gas have non-Gaussian statistics more pro- 
nounced at small scale. The extreme CVI (E-CVI) responsi- 
ble for the non-Gaussian tails of their probability density func- 
tions (pdf) form elongated coherent structures over 0.8 pc. These 
structures have been tentatively identified with regions of intense 
velocity- shear Q and enhanced local dissipation rate, based 
on their thermal, dynamical, and chemical properties. These 
pure velocity- structures do not follow those of dense gas, they 
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^ In the following, "shear" is used instead of gradient to empha- 
size that the observations provide cross-derivatives of the line-centroid- 
velocities (CV) (i.e. the displacement, in the plane-of-the-sky (POS), is 
perpendicular to the projection axis of the velocities). 



tend to be parallel to the magnetic field orientation, they are 
associat ed with gas warmer (Tv^^ > 25 K) than the bu lk of 
the gas (iHilv-Blant & Falgaronef[2007at iHilv-Blant et"ani2008l 
hereafter Paper I and Paper II), and they bear chemical sig- 
natur es ascribed to a warni chemistry not drive n by UV pho- 
tons (iFalgarone et al.1 120061: iGodard et aP l2009h . Last, in one 
such E-CVI structure, Plateau de Bure Interferometer (PdBI) 
observations disclose several sub- structures of intens e velocity- 
shear at scales as small as 6 milli-parsec (mpc) (Ealgarone et al.l 
I2009L hereafter EPH09). This suggests that turbulent molecular 
clouds harbour coherent velocity- shear structures from 6 mpc to 
800 mpc. 

We have increased the dynamic range by a factor 8 com- 
pared with Paper II, by mapping four times larger an area in the 
Polaris Elare with twice better a spatial resolution. The aim is 
to further explore the range of scales over which the spatial co- 
herence of these intense-shear structures is found. These are the 
first large scale observations performed at high-angular resolu- 
tion and high spectral resolution in translucent gas. The observa- 
tions and the results are described in Section 2 and 3. We briefly 
discuss possible interpretations and the nature of these structures 
in light of recent numerical simulations (Section 4). 

2. Observations 

Observations of the ^^C0(2 - 1) line were carried out at the 
IRAM- 30m tele scope with the 1.3m m multibeam heterodyne 
receiver HERA (ISchuster et al.ll2004 during August 2007 and 
January 2008, under good weather conditions. The map covers 
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R.A. Offset [arcsec] 

Fig. 1. Left: Integrated intensity map (Kkm s"\ T^scale) smoothed to 15". The white dashed boxes show the areas used to build the 
average p - v cuts (right panels). The blue rectangle (dashed line) delineates the previous field of Paper I. Right: p - v cuts in the 3 
boxes shown in left panel (r^ scale, distance measured in " with arbitrary origin). The line CV are shown in red. 



0.3 deg^ and consists of 9 submaps of 10' x 10', each of which 
was observed in two orthogonal scanning directions to minimize 
striping due to gain or atmosphere variations. T he final map en- 
compa sses the fields successively observed by iFalgarone et al.l 
(Il998h and in Paper I shown as a dashed box in Fig[T] and is 
~ 2 X 1 pc at the adopted distance of the source (150 pc). 
Data were acquired in the powerful on-the-fly (OTP) frequency- 
switched mode (4" scanning velocity, l^- time sampling, 4" 
spatial sampling in both directions, 13.8 MHz frequency throw), 
using the VESPA autocorrelator facility as backends. A total of 
1.5x10^ raw spectra was recorded in 80 hours of telescope time 
with a spectral resolution of 0.05 km s"^ Data were reduced with 
the ne w CLASS90 software optimized for OTP (iHily-Blant et al.l 
I2005h . The instrumental response was canceled by subtracting 
linear baselines to each original spectrum, which were then con- 
volved by a Gaussian kernel (1/3 HPBW) and gridded on a reg- 
ular grid with 0.5HPBW sampling. The final data cube was 
then smoothed to 15" and 0.1 kms"^ resolutions to improve 
the signal-to-noise ratio. The typical rms in each final pixel is 
Icr = 0.5 K in 0.1 km s"^ channels. 

3. Results 

3.1. Space and velocity maps 

The ^^C0(2 - 1) integrated emission is displayed in Pig. [T] (left 
panel) with three position- velocity (p - v) diagrams (right pan- 
els) made along the NE-SW boxes shown. A sharp variation of 
velocity, from ~ -3 kms"^ to ~ -5 kms"^ (from NE to SW) 
occurs over a layer thinner than a few 100" in projection. Pig El 
that displays the integrated emission in two adjacent velocity in- 
tervals, at high (HV) [-3.5, -0.5] kms"^ and low velocity (LV) 
[-6.5, -3.5] kms"\ stresses a remarkable characteristic of that 
field: the edges of the LV and HV components follow each other 
closely in projection over more than ~ 1 pc. It is then most un- 
likely that they be unrelated pieces of gas along the line of sight: 
they have to be in contact. 



The second remarkable characteristic of that field is the fol- 
lowing: while the emissions in the LV and HV components are 
extended, their spatial overlap (the pink areas of Pig El also ex- 
emplified in the p - v diagrams) is limited to narrow filamen- 
tary regions in projection. It is most visible between 6 = 87°45' 
and 88° (thus over ~ 1 pc) where it does not split into several 
substructures. Now, if the LV and HV components were parsec- 
scale volumes, their interface would appear thin over ~ 1 pc 
only if viewed edge-on (within ±5° for a projected size less than 
one tenth of its real size), a case that we rule out on statistical 
grounds. We therefore infer that the ^^C0(2 - 1) HV and LV 
components are layers rather than volumes and that their inter- 
face is 1 -dimensional rather than 2-dimensional. This ensures 
that under any viewing angle the two extended velocity compo- 
nents present a narrow interface. 

The slope of the variations of the line CV drawn on the p-v 
diagrams provides a measurement of the velocity shear between 
the two components. On each cut shown, there is an average 
shear of 13 km s"^ pc"^ (a velocity variation of 2 km s"^ over 
^0.15 pc). Steeper slopes are also visible and locally provide 
higher shears up to 30 kms"^pc"^ (1 kms"^ over 0.03 pc) in 
the middle cut. These values are more than one order of mag- 
nitude larger (within the uncertainties due to projections) than 
the average value of 1 kms~^ pc~^ estimated at th e parsec scale 
in molecular clouds (iGoldsmith & Arquillall 19851) . The velocity 
field therefore significantly departs, at small scales, from pre- 
dictions based on the generally adopted scaling laws in molec- 
ular clouds. If velocity fluctuations over a scale / increase as 
Svi oc /^/^, velocity- shears should increase by no more than 
33^/2 ^ 5.7 between 1 pc and 0.03 pc. 

A closer inspection of the - v diagrams shows that (i) 
the sharpest variations of the line CV occur between two line- 
wings appearance (above -2.0 km s"^ for the HV wing and be- 
low -5.5 km s"^ for the LV wing), ( ii) the separation between the 
LV and HV wings steepens from top (0. 1 pc) to bottom (0.03 pc), 
and ( Hi) the layer of largest velocity shear coincides with the lane 
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Fig. 2. ^^C0(2 - 1) integrated intensity in two adjacent velocity 
ranges: [-6.5:-3.5] kms"^ in blue and [-3.5:-0.5] kms"^ in red. 
At a distance of 150 pc, 20' correspond to 0.9 pc. 



of enhanced ^^C0(2 - 1) emission visible in Fig. [T] at the center 
of each cut. 



3.2. Distribution of E-CVI 



Following iLis et al.l (Il996h : iPetv & Falgarond (l2003h and 
Paper II, we have built the pdf (see Fig. O of ^^C0(2 - 1) 
CVI over the whole field. The statistics is significantly improved 
compared to previous works. The probability density in the most 
extreme bins reaches 10"^. Fig. [3] displays the locus of the E- 
CVI. It is remarkable that the thin structure delineated by the 
i2CO(2 - 1) E-CVI in the SE area (blue box of Fig. [B is so sim- 
ilar to that obtained with the same method applied to a much 
smaller sample observed in a diff'erent transition, ^^C0(1 - 0) 
(Paper II). The E-CVI structure is the high-angular resolution 
view of the structure obtained with the same statistical analy- 
sis performed on KOSMA maps of the field (HPBW=120'0 and 
shown in Fig. 12 of Paper II. 

The E-CVI structure does not follow the peaks of the 
^^C0(2 - 1) line integrated emission. Instead, it coincides, over 
the ~ 1 pc region discussed in Section 3.1, with the narrow inter- 
face of the HV and LV components i.e. the intense velocity- shear 
(Fig.O center), and follows in detail the thin elongated structure 
in the extreme velocity range [-2.0, -0.5] km s"\ that of the red 
line wings (Fig.O right). This association between CO line wings 
and intense velocity- shears extends the findings of Paper II to 
higher-resolution and over a larger scale. 

These properties of the locus of E-CVI unambiguously sup- 
port, for the first time, the proposition of Paper II that the 
E-CVI trace intense velocity- shears in turbulent gas and that 
the extreme variations of the line CV are driven by the ap- 
pearance/disappearance of linewings over small scales. It is 
also the first observational proof of the early conjecture of 

garone &Phillipsl(ll990h that the broad CO linewings trace 
the intermittency of turbulence in molecular clouds. These re- 
sults clarify, at least in the case of translucent clouds, the 
controversy on the origin of small-scale CV variations ex- 
pected to be due primarily to density fluctuations, line-of- 



sight projections, and radiative transfer ( Esquivel et al.l l2007l: 
iMiville-Deschenes et al.ll2003l:lLevriell2004tK 

Last, this E-CVI structure is coherent over ~ 2 pc while 
its thickness is as small as 0.03 pc. Its aspect ratio is there- 
fore ~ 70. Its length seems to be limited by the size of the 
map (see the longer structure computed from the KOSMA data 
in Paper II). We also note that the E-CVI structure splits into 
multiple branches in several areas, in particular around off'sets 
(-1000", 800") and (-700", 500"). 

4. Discussion 

4.1. What is the nature of the interface? 

The interface is primarily an intense velocity- shear. The p - v 
diagrams show that this velocity shear corresponds to a discon- 
tinuity in the CO flow: the HV (LV) component is not detected 
above ~ 0.5 K in the SW (NE) of the shear. The flows undetected 
in the ^^C0(2 - 1) line are either CO-poor and/or too dilute to 
excite the transition. In this framework, we observe the yield of a 
strain developing in a gas undetected in the ^^C0(2 - 1) line: the 
gas we detect (denser and/or richer in CO) is generated in the 1- 
dimensional intense-shear interface and is spread in the POS by 
motions whose velocity cannot be measured. This scenario nat- 
urally produces the two components of the large velocity- shear, 
with sharp edges closely following each other over ~ 1 pc and 
little overlap in projection under any viewing angle. 

The intense- shear structure may however belong to a shock 
of unknown velocity in the POS. We have searched for SiO(2 - 1) 
line emission as a chemical shock signature within this struc- 
ture and found no emission above a significant low threshold 
3(T = 5 mK that corresponds, in the optically thin case, to a 
tiny SiO column density of about 10^^ cm"^. Hence, there is no 
chemical signature of C- shocks faster than 20 kms"^ detected 
at the scale of 0.03 pc (iGusdorf et al.ll2008b . But we cannot rule 
out a weak C-shock component (v^ < 2 kms"^) in the POS 
that would produce the density enhancement and/or the CO en- 
richment of the gas required to explain the non-detection in the 
^^C0(2 - 1) line of the gas before it enters the shear layer. This 
would be consistent with the sub- to trans-Alfvenic nature of the 
turbulent motions in that field (see Paper II). The solenoidal con- 
tribution to the interface (2 kms"^) would exc eed the possible 
compressive one (< 2 km s"^), in agreement with lFederrath et all 
( 2009) who find that our observed statistical properties of tur- 
bulence in the Polaris Flare are in very good agreement with 
solenoidal forci ng at large scale. This result is reminiscent of 
the finding of iMininni et alJ (l2006ah that the stronger the 
shear at large-scale, the more intense the intermittency of 
velocity increments at small-scale. 



4.2. A plausible link with the dense cores 

The above findings are similar to those inferred from PdBI ob- 
servations (FPH09) of the small (1' by 2') field shown in Figl3] 
(left). Velocity- shears as intense as 500 km s"^ pc"^ are detected 
there over distances of 6 mpc, at the edge of CO structures of 
velocities diff'ering by sever al km s~^. The PdB I field is close to 
two low-mass dense cores (iHeithausenI |2002|) , interestingly lo- 
cated at the tip of the E-CVI structure (Fig. [2 left). 

The viscous dissipation rate of tu rbulence being proportiona l 
to the square of the rate-of- strain (iLandau & Lifschit3ll987l) . 
it is tempting to interpret the large increase of the velocity- 
shear, from 30 kms"^pc"^ (Fig. [T]) to 500 kms"^pc"^ in the 
PdBI field, as due to the development of an instability in the 
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Fig. 3. Left: Map of the CVI (top panel, color scale in kms"^) computed for a lag of 4 pixels or 60", and normalized pdf 
(circles, bottom panel) compared with a G aussian distribution (cr = 0.19 kms"^, red). The red crosses (top pane l) indicate 
the position of the dense cores from lHeithau sen ( 200^. The rectangle delineates the PdBI field of iFalgarone et al.l (I2QQ 9). Middle: 
E-CVI (blue contours) overplotted on the integrate d intensity of Fig. [T] Right: E-C VI (blue contours) overplotted on the intensity 
integrated in the red wing interval [-2:-0.5] km s"^ (iHily-Blant & Falgaronell2007ah . 



large-scale shear. The growth of the instability splits the shear 
into small-scale and more intense shears, thus increasing the lo- 
cal dissipation rate of turbulence by two-orders of magnitude. 
Clustering of small-scale structures of high strain-rate magni- 
tude (and therefore large dissipation) into structures of inertial 
extension have been found in numerica l simulations of incom- 
pressible HP (Moisv & JimenezI |2004|) and MHD turbulence 
(iMininni eTa l. 2006b). One then may speculate that these bursts 
of dissipation eventually lead to the formation of low mass dense 
cores largely devoid of turbulent energy, after an evolution of the 
gas that remains to be understood. 

5. Conclusions 

We have detected a 1 -dimensional structure of intense velocity- 
shear (~ 15 to 30 km s"^ pc"^), coherent over ~ 1 pc with a thick- 
ness of only 0.03 to 0.15 pc. This remarkable structure follows 
the distribution of extreme ^^C0(2 - 1) line- wings and coincides 
partly with the locus of E-CVIs in the field. These findings sup- 
port the previous claim we made that, in translucent molecular 
clouds, E-CVIs are tracers of extreme velocity- shears in inter- 
stellar turbulence, as do the broad CO linewings. 

This shear structure is proposed to be the source of layers of 
CO-rich dense gas in a CO-poor (and/or dilute) gas component 
experiencing the strain and not seen in ^^C0(2 - 1). The shear is 
likely to be the site of enhanced turbulent dissipation. We cannot 
rule out an undetected shock component in the POS. 

These results, in conjunction with the PdBI results of FPH09, 
stress the coupling of small and large scales in interstellar tur- 
bulence, over a dynamic range never reached before, i.e. from 
6 mpc to more than 1 pc. They support a framework in which 
trans-Alfvenic (but supersonic) turbulence dissipates primarily 
in intense- shear layers connecting the large- scales to mpc scales 
(or below). 

We speculate that turbulence dissipation has been proceeding 
for a longer time at the southern tip of the E-CVI structure than 



in the northern part, leading to the formation of the two dense 
cores. 
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